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Abstract 

In recent studies the truncated Levy process (TLP) has been shown to be very 
promising for the modeling of financial dynamics. In contrast to the Levy process, 
the TLP has finite moments and can account for both the previously observed excess 
kurtosis at short timescales, along with the slow convergence to Gaussian at longer 
timescales. I further test the truncated Levy paradigm using high frequency data from 
the Australian All Ordinaries share market index. I then consider, for the early Levy 
dominated regime, the issue of option hedging for two different hedging strategies that 
are in some sense optimal. These are compared with the usual delta hedging approach 
and found to differ significantly. I also derive the natural generalization of the Black- 
Scholes option pricing formula when the underlying security is modeled by a geometric 
TLP. This generalization would not be possible without the truncation. 
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1 Introduction 



It has been widely appreciated for some time that fluctuations in financial data show consis- 
tent excess kurtosis indicating the presence of large fluctuations not predicted by Gaussian 
models. With the continuing growth in the derivatives industry, and the recent emphasis on 
better risk management practices, the need for models that can describe these large events 
has never been greater. 

Several authors jl] have explored the stable Levy class of distributions as a possible 
alternative to the Gaussian. Levy distributions exhibit scaling or fractal properties which 
occur commonly in complex systems often studied in statistical physics. Indeed the Levy 
distribution does seem to provide a consistently better representation of financial data than 
the Gaussian distribution. Despite this it would be fair to say that there has been a lack 
of interest in this class of models, especially in the option theory literature (see for some 
exceptions). One reason for this would have to be a reluctance to accept the infinite variance 
which is typical in these models. This makes it difficult to find an appropriate generalization 
of the firmly entrenched Black-Scholes option pricing framework. Perhaps a more important 
reason is that financial data tends to become more Gaussian over longer timescales ||. |J. 
This property is also evident in the decay of the implied volatility smile obtained from 
observed option prices with increasing maturity ||. These properties cannot be explained 
by Levy distributions due to their stable additive property (the central limit theorem does 
not apply due to their infinite variance). 

A much more popular approach to explaining the excess kurtosis is based on the observa- 
tion that the variance (or volatility in financial language) of financial data appears to behave 
randomly. Stochastic volatility will generate kurtosis in an otherwise Gaussian process and 
this has led to a body of literature that attempts to model volatility as a diffusion process 
(for a review see ||). An important reason for the popularity of this approach is that it is 
still based on a Gaussian framework making it a relatively Black-Scholes friendly explana- 
tion for excess kurtosis. A major failing of stochastic volatility models is that they do not 
describe the ubiquitous power law or scaling properties observed in financial data. These 
properties are observed in the volatility correlation function ||, the PDF of high frequency 
price increments M, [7| , and the temporal decay of the peak of the PDF describing the finan- 
cial process 0. These last two properties are in fact, for short time horizons, well described 
by a simple Levy process. 

The growing empirical evidence of power law properties in financial data has generated 
renewed interest in the Levy paradigm. Recent work has shown that the problems associated 
with the Levy distribution can be simply overcome by what is known as a truncated Levy 
distribution (TLD). The TLD is Levy like in the central part of the distribution, but has a 
cutoff in the far tails that is faster than the Levy power law tails. The cutoff will ensure the 
variance of the TLD is finite. Financial prices over time can be described with the truncated 
Levy flight (TLF) or its continuous time limit the truncated Levy process (TLP). The TLF is 
constructed from sums of independent and identically distributed random variables described 
by a TLD. Since the TLD has finite variance, the central limit theorem applies and the TLF 
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approachs a Gaussian distribution as desired. However what is interesting is the existence 
of a characteristic timescale separating the Levy and Gaussian regimes. This timescale 
can be arbitrarily long due to the stable nature of the Levy distribution. Mantegna and 
Stanley || were the first to make the above observations regarding the TLF. They drew 
their conclusions based on the sums of independent Levy distributed random variables with 
a discrete cutoff in the tails. Inspired by these results, Koponen |J derived an analytical 
form for the characteristic function of a TLD with an exponential cutoff in the tails. These 
theoretical results have led to several recent empirical studies all supporting the TLF as a 
simple and effective model of financial data f|, |?J. 

The accurate modeling of financial price series is important for the pricing and hedging 
of financial derivatives such as options. Research on option theory with alternative pricing 
models has tended to focus on the pricing issue. It is now well known that non-Gaussian 
pricing models lead to the familiar volatility smile effect caused by the 'fat' tails of the non- 
Gaussian PDF's. These effects are well known will not be the focus here. What is much less 
understood and discussed is the issue of option hedging for non-Gaussian models. This is 
surprising because this would seem in many cases to be a more important issue than pricing. 
For liquid options (usually the vanilla types discussed here), the price will be determined by 
the market. The real use of the model is to define a hedging strategy. The models option price 
simply provides a way of testing the model against the market price. The standard approach 
to option pricing and hedging is the Black-Scholes framework. What is remarkable about 
the Black-Scholes case is that, for Gaussian or log-Gaussian pricing models, there exists a 
hedging strategy which will eliminate all the risk to an option seller. This leads to the well 
known delta hedging result which says that the hedge value is given by the derivative of the 
option price with respect to the current price of the underlying security. However for more 
general pricing models a riskless hedge does not exist and in these cases the Black-Scholes 
framework does not tell us how to proceed. This is especially problematic for times close to 
expiry where deviations from Gaussian are large. The standard approach is to simply apply 
the delta hedging procedure regardless of the pricing model used. However this approach is 
purely ad hoc as it has no clear theoretical basis. 

The inadequacies of the Black-Scholes framework led Bouchaud and Sornette to develop 



a simple and more general approach to option theory nty, 11]. Although in general a riskless 



hedge does not exist, it is possible to find an optimal hedging strategy that will minimize 
some appropriate measure of risk. An obvious choice for a risk measure is the variance of the 
wealth distribution of an option seller (or its 4th moment which would place more weight on 
the tails). Bouchaud and Sornette derived an expression for this optimal trading strategy 
which is valid for any non-Gaussian pricing model. In general the optimal trading strategy 
is not given by the delta hedge, though it is recovered for the special case of Gaussian or 
log-Gaussian models. 

The outline of this paper is as follows. In section 2 we will further test the TLP model 
using high frequency data from the Australian All Ordinaries share market index. With 
the derived parameter values for the TLP, we compare in section 3.1 the Bouchaud- Sornette 
optimal hedging strategy, the TLP delta hedging strategy and the Gaussian delta hedging 
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strategy. Although the TLP has finite variance, in the early Levy dominated regime a 
tail distribution based method to find an optimal trading strategy may be preferable to a 
moment based method. With this in mind, in section 3.2 we adapt to the TLP a simple 
tail distribution based method used by Bouchaud et-al 0, ^° ^ n( ^ the optimal hedging 
strategy for the Levy process. Finally in section 4 we will derive a natural generalization 
of the Black-Scholes option pricing formula for the case when the underlying security is 
modeled by a geometric TLP. This demonstrates that, unlike for the plain Levy process, the 
Black-Scholes framework is easily adapted to the truncated Levy paradigm. 



2 The Truncated Levy Distribution 

In this section we will outline the essential properties of the TLD and the TLF. Consider a 
general probability density function (PDF) P(x) and its characteristic function (CF) P(k) 
defined by 

1 fOO A A poo 

P( x ) = — dk P(k)e ikx , P(k) = / dx P(x)e~ ikx = (e- ikx ). (2.1) 

27T J —oo J —oo 



Moments of the distribution can be found from the CF by using 



/ flm N 

<* m W> = K^(M) 



fc=0 



(2.2) 



Normalization of the PDF requires P(0) = 1. We will deal purely with symmetric P(x) 
which in turn requires P(k) to be real and symmetric. 



The symmetric Levy distribution is defined by the CF [12 



L(k) = exp(—c a \k\ a ), 0<a<2 (2.3) 

where c is the scale factor and a is the characteristic exponent. The full PDF for the 
Levy distribution is only known analytically when a = 1 (Cauchy distribution) and a = 2 
(Gaussian distribution). However the value of the Levy distribution is known at the origin 
where 

L(x = 0) = ^ (2 .4) 

nac 

and in the tails where (for a < 2) 

c a r(l + a)shi7ra/2 

L{x) -> — — , x -> oo. 2.5 

7r|x| a+i 

These 'fat' power law tails mean the fractional moments (ja;^) are finite only for fi < a. In 
particular, for a < 2 the variance is infinite. 
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The TLD is a generic description for a Levy distribution that has some cutoff far in the 
power law tails. Such a cutoff will ensure that the variance of the distribution is finite. One 
possible cutoff is the exponential function for which the CF has been shown to be |J 



T(k) = exp 



k 2 + \ 2 ) a/2 cos {a arctan(£;/A)} - A c 



cos(7ra/2) 

With this CF the TLD in the tails takes the form 

c a T(l + a) sin(7ra/2)e- A l :E l 



(2.6) 



Tlx) 



1T\X 



l+a 



X 



OO. 



(2.7) 



Clearly the TLD reduces to the Levy distribution when A = 0. The asymmetric general- 
ization of (2.6) was also derived in ||. However in the paper we will deal only with the 
symmetric version. Applying (2.2) to (2.6) we find the variance and kurtosis of the TLD to 
be 



all — a) 



a 



a \ a— 2 



■c a \ 



k 



cos(7ra/2)(a: — 2) (a — 3) 



cos(7ra/2) 

where the variance and kurtosis are defined by 



a(l - a)c a X a 



a 



(x 2 ), k 



(x 2 ) 2 



3. 



(2.8) 



(2.9) 



We can always set c = 1 by scaling iasa;-> 'jx. The CF of the scaled x is again (2.6) but 
now with c — ► cy and A — ► A/7. The kurtosis (2.8) and the exponent a will both remain 
fixed under scaling. 

An important consequence for option pricing is that the exponential cutoff in (2.7) will 
ensure that exponential moments (e nx ) (n is any real number) exist for n < A. The expo- 
nential moments can be found from the CF (2.6) simply by substituting k 2 - 



-n 



This 

is equivalent to putting k — > in in (2.1) when the CF is symmetric and can be written as a 
function of k 2 . After the substitution we can expand (2.6) in powers of n 2 /X 2 and find to 
first order that 



(e nx ) ~ exp 



n 



{x 



A 2 » n 2 . 



(2.10) 



In the Gaussian case (A = 0, a = 2) this is exact. 



2.1 Convergence to Gaussian 

Let us now consider x to be the sum of iV independent and identically distributed random 
variables Xj with a TLD defined by the CF (2.6). Following Mantegna and Stanley || we 
can refer to a; as a TLF. The CF of x will be 



f(k, N) = exp 



c a N 



cos(7ra/2) 



{{k 2 + X 2 ) a/2 cos {a arctan(A;/A)} - A Q ) 



a^l (2.1i; 
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which is the same as (2.6) but with c a — > Nc a . First consider the special case when each Xi 
has a PDF described by a Levy distribution (A = 0). The Levy distribution is also known as 
the Levy stable distribution. This is because the PDF of the rescaled variable xN~ l ^ a is the 
same as that of Xi, that is the Levy distribution is stable under addition (the central limit 
theorem does not apply since the Levy distribution has infinite variance). We can also refer 
to this as the scaling or fractal property of the Levy distribution. On the other hand, for 
A > 0, the variance of Xi is finite so the PDF of x must approach a Gaussian by the central 
limit theorem. What is remarkable is that, as first pointed out by Mantegna and Stanley 
H, the convergence of x to a Gaussian occurs extremely slowly due to the stable property 
of the Levy distribution. 

We can derive the crossover time N c for x to converge to a Gaussian as follows. From 
(2.8) we know that the variance and kurtosis of x are given by 

a\N)= a{ )-°;\ j*\°-*N, fc(iV) = ^(W2)(«-2)(«-3) 

V ; cos(7ra/2) ' v ; a(l-a)c a X a N v ; 

A useful qualitative model of the distribution of x is that of a Gaussian in the central part 
out to the scale determined by the square root of the variance (2.12) (and therefore growing 
as \/~N). Beyond this scale we can think of the distribution being described by the tails (2.7) 
(with c a — > Nc a ) which are slowing being consumed by the Gaussian part. We can then ask 
at what time N c does the scale set by the square root of the variance equal the cutoff scale 
A -1 . We easily find the crossover time N c to be 

= c-JWa^ (2.13) 
a(l — a) 

A more rigorous method to derive this timescale is to find the time for the kurtosis of x to 
decay away. From (2.12) we see that this timescale is consistent with (2.13) which justifies 
the simple qualitative picture of the convergence to a Gaussian. We will find this picture 
helpful later when option hedging is considered. Mantegna and Stanley || first derived 
the timescale N c ~ c~ a X~ a using a method based on the probability of x returning to the 
origin (this method gives a slightly different a dependent coefficient). This result can also 
be obtained from the Berry-Esseen theorem Jl3| and will be independent of the precise form 
of the cutoff. Clearly N c can be as large as one wishes by making A small enough. 



2.2 Parameter Fitting 

Parameter fitting to a TLF is no more complex than fitting data to a Levy distribution. 
The parameters a and c can be obtained by fitting a Levy distribution to price increments 
(or log increments) on the smallest timescale such that linear correlations are negligible. A 
good rule of thumb is 30 minute increments. The timescale needs to be small so that we are 
operating well in the Levy regime of the TLF. Once c and a are known we can extract the 
cutoff parameter A from the variance of the dataset using (2.12). This will fit the variance 
of the data for a particular time N. From the point of view of option pricing fitting to the 
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daily variance would seem the most appropriate choice. It would be appropriate to fit A 
using a large dataset to ensure a good description of the tails. Once a and A are determined 
we can expect these parameters to be relatively stable despite the well documented non- 
stationary behavior associated with a stochastic variance. These effects can be accounted 
for by a time dependent scale factor c. The non-stationary behavior is associated with 
variance measurements over small data sets (typically 1 month). This data will only sample 
the central part of the PDF which compared to c is relatively insensitive to A. The relative 
stability of the exponent a over time might be expected since it is invariant to scale changes. 
As a consequence sums of Levy distributed random variables with the same exponent a, 
but different scale factor c, will still have a distribution described by the exponent a. This 
stability has been confirmed in empirical studies on the S&P 500 index [|J. 

There are several methods for fitting data to a Levy distribution (see Rachev and Mittnik 
for a review). Here we adopt the method used by Mantegna and Stanley ||. It is based 
on the relation 

T(x = 0, N) ~ L(x = 0,N) = ^{$ a , iV«iV c (2.14) 

which will be a good approximation in the Levy regime of the TLF (N <C N c ). In this 
method we can extract c and a by fitting a power law decay to the probability of return to 
the origin for times N <C N c . Alternatively we can take the log of (2.14) and extract c and 
a by a straight line fit. This is a very practical method because it is simple, it can be used 
on relatively small datasets (since T(0, N) is the maximum of the distribution), and because 
of the easy availability of high frequency financial data. It is somewhat different to other 
methods in that it fits the peak of the distribution at various times rather than fitting the 
whole distribution at a single time. 

We will demonstrate the parameter fitting process using high frequency equity data. We 
have used the All Ordinaries Index (AOI) which is the major Australian equity index. The 
dataset comprised the value of the AOI at 5 minute intervals from 1993-May/1997. From 
this, datasets describing the raw change in the AOI were constructed for 30 minutes, 1 hour, 
2 hours, 3 hours and 1 day. The mean was then subtracted from these datasets. The time of 
30 minutes was considered to be a good minimum time for which changes in the AOI could 
be considered to have negligible linear correlations. This is necessary in order to justify the 
use of (2.14). From these datasets the probability of zero change in the AOI was found and 
plotted against time in figure 1 (time is in units of 30 minutes with 1 trading day equal to 
6 hours or N = 12). 

The data of figure 1 was very well fitted to a power law decay curve. From the curve of 
best fit and (2.14) it was found that a ~ 1.2 and c ~ 1.1. For comparison in figure 1 the fit 
of -P(O) to a Gaussian process is shown where the Gaussian is fitted to the daily variance of 
226 (as would occur in practice from an option pricing perspective) and extrapolated back 
to 30 minutes. In figure 2 the empirical PDF for the 30 minute data is shown along with the 
Levy PDF defined by the CF (2.11) with the parameters a = 1.2, c = 1.1, A = and N = 1. 
We also show the Gaussian PDF which has been fitted to the variance of the daily data and 
extrapolated back to 30 minutes. The Levy PDF can be seen to be in good aggreement with 
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the empirical PDF all the way out to to ±50 points which is approximately 12 standard 
deviations. The tails of the distribution are quite sensitive to the exponent a which we 
demonstrate by also showing the Levy PDF for a = 1.4. From this we can see that simple 
fitting method used has been quite effective. These results along with others [[[], |4|, |7j further 
support the LF/TLF model of financial data. 

To fit A the dataset used was that of daily changes in the AOI between 1987-96. A large 
data set is required to fit A since this relates to the tails of the distribution. Substituting 
the daily variance of 345 and N = 12 into (2.12), and using the same values for a and c 
derived previously, we find that A = 1/80. With these parameter values we find that the 
Levy-Gaussian crossover time (2.13) is N c ~ 222, or approximately 19 trading days. This 
is very close to 1 trading month (21 days) which previous studies have shown to be the 
timescale for which the Gaussian description becomes accurate |3|, £|] . Since N c 3> 1 day the 
use of (2.14) out to 1 trading day is justified since this is well within the Levy regime of the 
TLF. 

The daily dataset from 1987-96 is quite symmetric except in the far tails. The largest 
positive deviation in the dataset is 80 points. On the other hand there are several nega- 
tive deviations beyond 100 points which all occured around the 1987 crash of -520 points. 
These facts suggest that an appropriate asymmetric model would be a symmetric Levy dis- 
tribution with asymmetric cutoff parameters, rather than an asymmetric Levy distribution 
with symmetric cutoff parameter. The later model || would give asymmetry in the whole 
distribution rather than just in the far tails. In option pricing generally only part of the 
probability distribution is important. This suggests that a simple and effective approach to 
asymmetry would be to choose a cutoff parameter that best fits the part of the distribution 
that is relevant. To illustrate this we have derived the cutoff parameter A + = 1/41 from the 
variance of positive deviations (203), and the cutoff parameter A_ = 1/122 derived from the 
variance of negative deviations (487). This latter cutoff parameter leads to a crossover time 
of approximately 30 days. We have used a = 1.2, c = 1.1 and N = 12 as before. 

3 Optimal Option Hedging 

In this section we will be concerned with the optimal hedging of call options for times to 
expiry which are less than the Levy-Gaussian crossover time. With 1 trading month as a 
typical crossover timescale we would be interested in times to expiry roughly less than 10 
trading days. In this case the drift and probability of negative prices are negligible and we 
can consider the TLP as our model for the financial data. 

Consider first the price C(S , E, t) of a European call option at current time t = 0, with 
exercise price E due to expire in a time t. When the time to expiry is small the returns and 
interest rates can be neglected. The option price is then very well approximated by 

C(S , E, t) ~ (max(S - E, 0)) = J™ dS (S - E)P(S, t\S , 0) (3.1) 

where P(S, t\So, 0) is the driftless PDF of the underlying asset price S(t). We know in 
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advance that (3.1) will lead to the familiar volatility smile effect caused by the 'fat' tails of 
non-Gaussian PDF's like the TLD. These effects are well known and we will not consider 
them further. What is much less understood and discussed is the issue of option hedging for 
non-Gaussian models. It is to this that we now turn. 



3.1 Variance based hedging strategy 

A framework which can address the hedging issue is the Bouchaud-Sornette approach to 



option pricing and hedging flnj, |11|] . This approach starts by finding the variation in wealth 



for a call option seller. They find the wealth variation between times and t can be written 

as 

AW\l = C(S ,E,t) -max(S(t) - E,0) + [ dr (j)(S T )S(r) (3.2) 

Jo 

where the first term is the option premium received at t = 0, the second term describes the 
payoff at expiry t and the third term describes the effect of trading where <f>(S T ) (S T = S(r)) 
is the amount of stock held. The interest rate is set at zero for clarity which will be a 
good approximation for short-medium term options. The option price is found by requiring 
(AW\ ) = 0. This leads to 

C(S Q , E, t) = (max(S(t) - E, 0)) - J* dr (<f>(S T ))(S(r)) (3.3) 

where the increment S(r) is in the future and assumed to be independent of S(r) (this 
assumption can be relaxed [fLOfl ). What is remarkable about the Black-Scholes case (Gaussian 
or log-Gaussian models) is that there exists a trading strategy 4>*{S T ) such that AVT|q can 
be made to vanish. In this case the rate of return dependence in (3.3) cancels and the Black- 
Scholes result follows [|K| . For more general pricing models a riskless hedge does not exist and 



in these cases the Black-Scholes framework fails. However the Bouchaud-Sornette method 
is easily adapted to these more realistic situations. Although in general a riskless hedge 
does not exist, it is possible to find an optimal trading strategy <f)*(S T ) that will minimize 
some appropriate measure of risk. An obvious choice for a risk measure is the variance of 
the wealth distribution (Aiy 2 |[)[0]), or its 4th moment which would place more weight on 
the tails. The optimal trading strategy can be easily computed for any PDF describing the 
price S(t). In these more general cases the option price (3.3) (evaluated at 0*) will depend 
on the rate of return. However the return can be safely set to zero for options less than a 
few months to expiry (see Aurell et-al [[(J). In this case the second term on the right hand 
side of (3.3) will vanish and (3.1) will be a good approximation. However the option price 
(3.3) will need to be corrected by a risk premium whose scale will be set by the residual risk 



(AH^Iol^*]). This risk premium can account for the bid-ask spreads in option prices. 
Consider the special case of a pricing model in which the price increments are uncorre- 
cted, stationary and have zero mean (these assumptions can be relaxed). This includes the 
TLP and will be a good model for options roughly less than a month to expiry. In this case 
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the hedging strategy that minimizes the variance of the wealth distribution has been shown 



to be PL 11 



^(S ,E,t) = ^- j°° dS (S-E)(S-S )P(S-S ,t), (3.4) 
a z t Je 

where P(S — Sq, t) is the conditional PDF of the pricing model with initial price So and a 2 t is 
the variance of the model. The optimal hedge 0* ranges between and 1 since we are dealing 
with an option on one unit of the underlying security. For at-the-money options (E = So) 
we easily find from (3.4) that <fi* = 1/2 for any pricing model. For options out-of-the-money 
(E > So) we have <fi* < 1/2 and for options in-the-money (E < So) we have <fi* > 1/2. 
For the special case of the Gaussian pricing model with PDF 

G(s - s °' tH 7fbr xp pi^)' (3 - 5) 

the optimal hedging strategy (3.4) reduces to the Gaussian Black-Scholes result 

<P*(So,E,t) = —J E dS (S - E)G(S - S ,t) = N l-j-^-j (3.6) 

where N is defined in (4.17). This result can also be shown to be the probability of exercise 
of the option. This hedging strategy is clearly just the derivative of the option price with 
respect to S . This is the well known delta hedging result and can be expressed more 
generally as 

A(5 , E, t) = -£rC(S , E,t) ~ ^- f° dS (S - E)P(S - S , t) (3.7) 
obo obo Je 

where we have used the approximate option price (3.1). This is a riskless hedging strategy 
only for Gaussian and log-Gaussian models. However it can be used to give the hedge value 
for other pricing models. 

Consider the PDF for the TLP which has the tail form 

rpfc c Q tr(l + a)sinvra/2 \S - S \ , . 

T{S-S ,t)~ 7r|g-S |i+° exp(-A|S- S \), -j= > 1. (3.8) 

The scale where this approximation becomes valid is set by the Levy process diffusion scale 
ct 1/Q . The variance of the TLP will always be greater than this scale for times less than the 
Levy-Gaussian crossover time. We will always be interested in times less than this which 
makes the condition on \S — So\ appropriate. We wish to compare for this model, the optimal 
Bouchaud-Sornette hedging strategy (3.4) with the delta hedge (3.7). We will also compare 
the TLP hedges with the Gaussian delta hedge (3.6). Applying the PDF (3.8) to (3.4) we 
find that the optimal Bouchaud-Sornette hedge becomes 

<KS 0t E, t) * + a) ^ \y 2 ~ a e- y + (1 ~ V - a)r(2 - a, y)} , > 1, (3.9) 



2ira(a — 1) : 
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where 



and 



/y=lrV /2! ^# 1 = A|^-5 | 



rfa, x) 



'<r 2 t 



(3.10) 



(3.11) 



is the incomplete gamma function. In (3.10), t c is the same as the discrete time Levy- 
Gaussian crossover time defined in (2.13). For So > E we have 



•(S ,E,t)~l 



T(l + a) sin7ra 
2na(a - l) 2 



y 2-a e -y + {l _ y _ a)r{2 _^ y) 



Sn-E 



> 1. 

(3.12) 

Consider next the delta hedging strategy based on the TLP. Applying the PDF (3.8) to 
(3.7) we find 



A(S ,E,t)~ 



t T(l + a) sin7ro; 
T c 2ira 2 (a-l) 2 



y l - a e~ y -(a- l)y- a e~ y - T(2 - a, y) 



E~S 



>1, 



and for S > E 



A(S ,E,t)~l- 



(3.13) 



t T(l + a) sin7ro; 
T c 27ra 2 (a-l) 2 



y L-° e -y _ ( a _ i) y 



-« e -w-r(2-a,s/) 



So — E 



> 1. 



(3.14) 

The Gaussian delta hedge (3.6) is exactly equal to the probability of exercise of the option 
for the Gaussian pricing model (3.5). The TLP delta hedge (3.13-14) is also equal to the 
probability of exercise when calculated with the tail approximation (3.8). 

In figures 3 and 4 we plot and compare the TLP Bouchaud-Sornette optimal hedge (3.9), 
the TLP delta hedge (3.13) and the Gaussian delta hedge (3.6). We are using E > So 
which means the options are out-of-the-money and have a low hedge value and probability 
of exercise. We plot the hedge value against the ratio of E — S and the standard deviation 
of the pricing model. No parameter values need to be specified in the Gaussian delta hedge 
(3.6). For the hedges (3.9) and (3.13), we know from (3.10) that we need to specify the 
exponent a and the ratio of the time to expiry t and the Levy-Gaussian crossover time t c . 
In section 2 we obtained a — 1.2 and a Levy-Gaussian crossover time of approximately 
20 trading days. We have used these parameters with the expiry times of 1 and 5 days 
in figures 3 and 4 respectively. In figure 3 we see that Bouchaud-Sornette optimal hedge 
leads to hedge values that are, in percentage terms, significantly greater than those obtained 
with both delta hedges. The relative difference increases as E — So grows. In figure 5 the 
differences are not as great, but in percentage terms they are still very large. 



3.2 Distribution based hedging strategy 

The Bouchaud-Sornette method will fail for models such as the Levy process which have 
infinite moments. Although the TLP has finite variance, in the early Levy dominated regime 
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a distribution based method to find an optimal trading strategy may be preferable to a 
moment based method. Below we will adapt to a truncated Levy model a simple method 



used by Bouchaud et-al 0, 11] to find the optimal hedging strategy for Levy models. 



Consider the case where <f> is the final hedge value some time t from expiry. Then from 
(3.2) we find that the change in wealth over this period is (ignoring the premium) 

AW = J((0 - l)(S t -S )+E- S ) + (1 - I)(f)(S t -S ), < <j> < 1 (3.15) 

where 1 = 1 for S t > E and 1 = for St < E. Now we make the key assumption that / can 
be treated as an independent random variable that takes on values 1 and with probabilities 

V and 1 — V, where V is the probability the option will be exercised. Note that both <fi and 

V are functions of (So,E,t). We can write the change in wealth AW\s t -s , due only to the 
change in price as 

AW\ St s = 1(1 - <P)(S - S t ) + (1 - I)(j){S t - So). (3.16) 

In most realizations of / and St — So, this will lead to a loss which we denote by I. We 
assume that St — So is described by a TLP which has the PDF (3.8) in the tails. Due to 
the independence assumption of I we can easily write down the tail PDF of the loss I using 
(3.8). From this we find that the probability of a loss greater than Z* is 



CO 



e -\-l/4> ,00 g-A + J/(l-0) 



v(i>u)~{i-v)r j h di^j^+ni-cpr j u * ll+a . (s.i?) 

We can write /* = 1/(2A*) which is the average loss suffered to an event ±1 / A*. Thus we can 
consider minimizing losses above Z* to be equivalent to minimizing losses to events beyond 
±1/A*. So after performing the integrals we can write (3.17) as 

V{1 > Q ~ (l - V(S , t)) \°LT(-a, p.) + V(S , t)X a + T(-a, p+) (3.18) 

where 

"- = 2^ = (3 ' 19) 

We wish to find the optimal trading strategy 0* that will minimize the probability (3.18). 
This is obtained by by taking the derivative of (3.18) with respect to <ft, setting the left hand 
side to zero to obtain 

0= (l-p)0 Q - 1 exp(-/3_)-P(l-0) Q - 1 exp(-/3 + ), (3.20) 

and solving this equation for <fi* as a function of V . For A = 0, (3.20) becomes independent 
of A* and it can be solved to obtain the Levy optimal hedge 

ffl^' ^wSw ' (3 - 21) 
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first derived by Bouchaud et-al |2|, [11]]. For a = 2 the Gaussian delta hedge 0* = V is 
recovered. 

The Levy hedge (3.21) is of fundamental interest as it shows that a hedging strategy can 
be derived even for a model with infinite variance. However what is required in practice is a 
trading strategy that interpolates between the Levy hedge and the Gaussian hedge over the 
crossover timescale (2.13). This follows from the central limit theorem which ensures that 
the optimal hedge tends to the Gaussian hedge (this is not the case for the geometric TLP 
discussed in section 4). Clearly we want to choose A" 1 to be the smallest possible value such 
that the tail form of the PDF will hold. The scale beyond which this holds is set by the 
Levy diffusion scale ct 1 ^. A good choice is A7 1 = ay/t which is the standard deviation of 
St — Sq. This scale will always be greater than the Levy diffusion scale for times less than 
the crossover timescale. This is fine since this method of hedging will only be appropriate 
for these times. The choice for A* is also consistent with the qualitative picture, discussed 
in section 2.1, for the convergence of the TLP to a Gaussian. 

In figure 5 we plot the optimal hedging strategy against the probability of exercise ob- 
tained by a numerical solution of (3.20). We have chosen as parameters those obtained 
in section 2.2 (a = 1.2 and A = 1/80) which gave a Levy-Gaussian crossover time of ap- 
proximately 19 trading days. Using A~ x = o\ft, we plot the optimal hedge at 2 days 
(A/A* ~ 0.33), 10 days (A/A* ~ 0.73) and 20 days (A/A* ~ 1.03) along with the Gaussian 
hedge [a = 2, A = 0) and the Levy hedge (A = 0). The hedging strategies have the very 
nice feature of a well defined evolution from the Levy hedge (3.21) to the Gaussian hedge 
in about 10 days. Between 10 and 20 days there is a much slower creep away from the 
Gaussian hedge. From this we can tentatively conclude that the hedging strategy defined by 
(3.20) could be reliably used out to time t* defined by A* ~ A/0.73, or approximately 1/2 
the crossover time (2.13). Beyond this the Bouchaud-Sornette variance based method would 
be appropriate. In figure 5 we see that for options with a small probability of exericise, 
the difference between the optimal hedge and the Gaussian Black-Scholes hedge is large in 
percentage terms. This is consistent with what was found in section 3.1. 

Figure 6 is similar to figure 5 except that we have used the asymmetric cutoff parameters 
A + = 1/41 and A_ = 1/122 discussed in section 2.2. The symmetric cutoff parameter 
A = 1/80 was used in the variance required to describe the growth in A*. This case gives an 
optimal hedge that involves holding less stock than in the symmetric case. This is what we 
expect since the asymmetric case has more weight in the negative tails and negative price 
changes favor holding zero stock. Plots of the ratio of the asymmetric/symmetric hedging 
strategies of figures 5 and 6 approximately range between 0.7 and 1 for exercise probabilities 
between and 1 respectively. 



4 Option Pricing 

Here we will derive a generalization of the Black-Scholes option pricing formula for the case 
where the underlying security is modeled by a geometric TLP. We will see that the Black- 
Scholes framework is easily adapted to the truncated Levy paradigm, a feature not shared 
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by the plain Levy or geometric Levy process. 

The Black-Scholes option pricing theory |L4| is based on the geometric Brownian Motion 
(GBM) paradigm of financial market dynamics 
define the geometric TLP S(t) by 



15]. A natural generalization of GBM is to 



x(t) + fit — 



S(t) = S exp 
where x(t) is a stochastic process defined by the CF 



f g (k,t) = exp 



cos(7ro;/2) 



{k 2 + A 



2W2 



cos {a arctan(A;/A a; )} — A") 



(4.1) 



a^l, (4.2) 



a 2 t is the variance of x(t) (a 2 is defined in (2.8)) and \x will be shown to be the rate of return. 
A time dependent variance can be considered by replacing a 2 t with J ds o 2 (s)ds. This time 
dependence will derive from a time dependent scale factor c x (t) in (4.2) with a and X x kept 
constant. GBM is recovered by setting a = 2 and X x = in the CF (4.2). Using (2.8) this 
gives 2c 2 , = a 2 . 

Using (2.10) we find the moments of the geometric TLP are 



(S n {t)) ~ exp nfjt + n(n - l)a 2 t/2 



A . » n 2 



(4.3) 



From this we see that \x is the rate of return. This expression is an approximation to 
first order in n 2 /X 2 and is exact for GBM. The approximation breaks down for very high 
order moments which is direct consequence of the absence of a central limit theorem for 
multiplicative processes like (4.1). But how do we know that the condition A 2 3> n 2 will 
hold? Consider short timescales where x(t) is small and we can write (4.1) as 



S(t) ~ S + S x(t). 



(4.4) 



Under this approximation S(t) reduces to the plain TLP. Equation (4.4) will be a valid 
approximation when V a 2 t <C 1. This will apply for time periods typically of the order of 1 
month or less. Fitting data to raw changes in the price S(t), as in section 2.2, effectively fits 
Sox(t) to a TLP which is described by the parameters a, c and A. We can then extract the 
parameters c x , X x of the CF of x(t), by the relations c x = c/ So and X x = So A (a is invariant 
under scaling). We know that X x must obey X x ^> 1 because the original cutoff scale A -1 of 
the TLP must be much smaller than So to ensure the price stays positive. We have found 
that fitting daily data to change in log price leads to X x typically in the range of 20-30. 

Lets now consider the price C(Sq, E,t) of a European call option at time t — with 
exercise price E due to expire in time t. We will assume that the price of the underlying 
asset S(t) follows a geometric TLP (4.1) with S the current price. A simple but ad hoc 
approach to option pricing is to apply the the risk-neutral approach ]15[ in exactly the same 
way it is used with the GBM model. In this case the option price if given by 



C(S , t) = e- rt (max(S -E,0)) = e~ rt J°° dS {S - E)T g {S, t\S , 0) 



(4.5) 
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where T^(S,t\S ,0) is the PDF of the geometric TLP with ji set equal to r. Using 
Jo° — Jo ' the normalization of T g (S, t\S , 0) and (4.3) we find that 



C(S ,t)~S -Ee- rt -e- rt [ dS (S - E)T g (S,t\S ,0), A 2 . > 1. (4.6) 

v 

We can write the PDF of S(t) in the form 

dx r°° 

T g {S,t\S ,0)dS = — J dkT g {k,t) cos kx, (4.7) 

where T g (k, t) is the CF of x(t) defined by (4.2), and from (4.1) we have (with the risk-neutral 
measure) 

x = \n(S/S ) -rt+ \oH. (4.8) 

We can substitute (4.7) into (4.6) and change the order of integration. Using the integrals 

f Xe , , sm kx e . ,, . , A n . 

/ dx coskx = — \-7io(k) (4.9) 

J — OO rC 

and 

f Xe j x i x (coskx e + A; sin kx e ) 

/ dx e x coskx = e Xe - — — (4.10) 

J —oo 1 ~\~ k 

we find the call option price is given by 

^v/ „ ^ x ^ 1 r-t Ee~ rt f°° - ,, . /sin /cx e — cos /cx e \ . 9 „ . , _ „ N 
C(S , J M^ --£e- rt + ^^y o dfcT,(M)^ ^ + fc2) "J , ^»1 (4.11) 

with x e defined by 

x e = \n(E/S ) - rt + i<r 2 t. (4.12) 

We found A x to be in the range 20-30. This means the approximation (4.11) will be very 
good. 

The Black-Scholes result is recovered from (4.11) by writing 

f g (k, t) = exp(-a 2 tk 2 /2) (4.13) 

and using the identities 

sin kx e — k cos kx e sin kx e (k sin kx e + cos kx e ) 



k{\ + k 2 ) k 1 + k 2 

x — 2 / ~fp^ 



(4.14) 



r°° e 9 fx — 2^ /, d\ 

/ dk— -(ksinkx + ^coskx) = ire p ' f ~ x m I , 7 >0 (4.15) 

Jo Y + « V v 2p / 
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and 



where 



nk 2smkx 7r 
dk e~ pk — - — = iiN 



2p, 



N(x) 



'2tt J- 



dye~ y2/2 , N(x) = 1 — N(—x). 



We then obtain the Black-Scholes option pricing formula |I3 



C(S , E, t) = S N(d) - Ee~ rt N(d - VaH) 



where 



d 



at — x t 



(4.16) 
(4.17) 

(4.18) 
(4.19) 



Its important to emphasize that the option price (4.11) will not approach the Black-Scholes 
result (4.18) as the time to expiry increases. This is because the central limit theorem does 



not apply to multiplicative processes like the geometric TLP [18 



5 Conclusion 

In this paper we have further demonstrated that the TLP can effectively model the em- 
pirical successes of both the Levy and Gaussian distributions at short and long timescales 
respectively. In option theory a major disincentive for using non-Gaussian based models is 
the absence of a riskless hedge. This makes it impossible to apply the Black-Scholes option 
pricing framework in anything other than an ad hoc way. We therefore discussed how the 
Bouchaud-Sornette approach to option theory provides a relatively simple solution to the 
problem of option pricing and optimal hedging in non-Gaussian models. We applied this 
optimal hedging strategy to the TLP pricing model and compared it with the delta hedging 
strategy for both the TLP and Gaussian models. Significant differences were found. We 
then discussed an alternative and computationally simple tail distribution based hedging 
strategy appropriate for the Levy regime of the TLP. For a TLP with exponential cutoff, the 
exponential moments exist and this allowed us to derive in section 4 a natural generalization 
of the Black-Scholes option pricing formula for the case of a geometric TLP pricing model. 
This is not possible for the geometric Levy process, an important reason for the lack of 
interest in this model. 

Further theoretical work could involve deriving CF's of TLD's with alternative cutoffs 
such as stretched exponential tails, or models with asymmetric cutoff parameters. Stretched 
exponential tails are interesting technically because the Laplace transform of the distribution 



would then exist [|T^]. An asymmetric cutoff will give a PDF that is symmetric in the center 
but asymmetric in the far tails. This was how asymmetry occurred in the daily data set of 
the AOI considered in this paper. Also of great practical importance is the construction of 
efficient algorithms for the numerical simulation of random variables with a TLD. It may be 
possible to adapt existing methods for the simulation of Levy distributed random variables 
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TTfl . With simulations we could test the effectiveness of the tail based hedging strategy 



against the variance based method. This is important since uncontrolled approximations 
had to be made in deriving the tail based strategy. The simple TLP model studied here 
does not address the apparent stochastic dynamics of the volatility along with its long range 
power law correlation function || . A challenge for the future is to construct pragmatic models 
which can describe these properties as well as the those of the distribution considered here. 
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Figure 1: Fit of Eq. (2.14) to high frequency data 
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Figure 2: Fit of the Levy PDF to 30 minute data using parameters derived from figure 1 
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Hedging Strategies: 1 day to expiry 
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Figure 3: Plot of the the TLP Bouchaud-Sornette optimal hedge (3.9), the TLP delta hedge 
(3.13) and the Gaussian delta hedge (3.6). 
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Hedging Strategies: 5 days to expiry 
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Figure 4: Same as figure 3 except 5 days to expiry 
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